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ABSTRACT 

We describe progress on a re-evaluation of the spectrum of cosmic 
rays determined with the Haverah Park shower array. Particular 
attention is paid to the reality of some giant showers. 

1, Introduction . We are engaged in a re-appraisal of the energy 
spectrum of cosmic rays above 10 18 eV as determined with the Haverah 
Park shower array. Here we offer a progress report on work which is 
motivated by the continuing controversy over the shape of the spectrum 
above 10 19 eV - in particular the Yakutsk group have questioned the 
reality of events of 10 2 ° eV - and by the recent re—inves tigation of the 
predicted shape of the spectrum above 10 19 eV if the sources of these 
particles are at cosmological distances (Hill and Schramm 1985) . 

2. Cosmic Ray Energy Spectrum . The differential energy spectrum 
derived from our work above 10 18 eV is shown in Figure 1. Above 
3.5xl0 18 eV the spectrum has been updated by the addition of events 
recorded to December 1983. A detailed analysis of possible sources of 
systematic error has been made taking into account the effects of un- 
certainties in zenith angle, lateral distribution fluctuations, core 
location and attenuation length, (Cunningham 1982). For energies 
between 8xl0 17 and B.SxlO^eV systematic selection effects and 
analysis errors dominate over statistical uncertainties and detailed 
simulations have allowed a deconvoluted spectrum to be derived. Above 
3. 5 x 10 1 8 eV the error analysis has been conducted on a shower-by-shower 
basis and the statistical errors have been shown to be at least twice as 
great as the instrumental errors. We do not yet regard the spectrum of 
Figure 1 as our ’final’ spectrum as further refinements will be possible 
as our detailed knowledge of showers increases but we wish to emphasise 
that we have considerable confidence in the durability of the intensities 
and energies assigned above 10 19 eV, The major differences between this 
spectrum and those published at Kyoto are (a) the exclusion of events 
with 0>45° (as we now regard our knowledge of the structure function to 
be incomplete above this angle) and (b) use of an energy dependent 
structure function measured in showers of 10 17 -5xl0 1 eV (Coy et al 
1981) and in a small number of large showers which fell during the period 
of that experiment. The main features of the spectrum are the flattening 
above 10 19 eV and its continuity to just beyond 10 2 ° eV . At the Paris 
conference we pointed out that the flattening may also be interpreted as 
a dip in the spectrum (Bower et al 1981) and suggested that if particles 
above a few times 10 18 eV were extragalactic origin then the dip might 
well be due to electron-pair production. This interpretation has been 
confirmed by the detailed analysis of Hill and Schramm (1985). 

The Haverah Park and Yakutsk spectra (Vaselev et al 1983) are compared 
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in the lower part of Figure 1. The spectra are found to agree reasonably 
well until about 3xl0 1 ^eV when the absence of large showers in the 
Yakutsk spectrum becomes apparent. 


In recent years a 
number of 1 m 2 blocks 
of scintillator have 
been incorporated in 
the Haverah Park array 
making possible a 
comparison between the 
model calculation con- 
version used by the 
Haverah Park group 
and the calorimetric 
approach of the 
Yakutsk group. We 
have shown else- 
where that the cali- 
bration is good (to 
within better than 
20%) up to at least 
5 x 10 19 eV and 
similarly that the 
Volcano Ranch energy 
estimates are in 
accord (Bower et al 
1983a, b). Also the 
Sydney experiment 
offers evidence of a 
flattened spectrum 
above ^ 4 x 10 19 eV 
(Horton et al 1983) ; 
that spectrum may 
extend to 4 x 10 2 ° eV 
(Linsley 1983). We 
do not plot the 
Sydney spectrum here 
because of un- 
certainties about the 
energy calibration. ^ 



3. Events of IQ 20 eV . The 4 most energetic events included in the 
spectrum have been assigned energies > 10 2 ° eV. Brief details of these 
are given in Table 1; maps of the density pattern observed in each event 
were published in the World Data Catalogue although the sizes have been 
slightly altered as a result of the revised lateral distribution function 
now adopted. Three of the events have risetime information available at 
one or more of the 34m 2 detectors and are discussed in that context in 
HE 4.7-6 (Lawrence et al) . 

Of the events in Table 1 by far the most outstanding in terms of number 
of densities and precision of core position is 17684312. Unfortunately 
this event was recorded in the epoch before scintillator densities were 
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being recorded. One of the most energetic events with scintillator 
density information is 21220296, a map for which has been published else- 
where (Bower et al 1983c) and these two events are contrasted in Table 2. 


Table 1 


Reference 

number 

angle 0 

a 

6 

b 

r } (m) 

Energy 

(eV) 

World Data 
Catalogue 

Rise- 

times 

8185175 

35 

353° 

19° 

1 

-O 

O 

o 

443 

1.02 x 10 2 ° 

p78 

None 

17684312 

35 

201° 

71° 

46° 

376 

1.05 x 10 2 ° 

p86,87 

1 

9160073 

30 

199° 

44° 

73° 

1384 

1.05 x 10 2 ° 

p79 

2 

12701723 

29 

179° 

27° 

o 

00 

1093 

1.21 x 10 2 ° 

p83 

4 


Table 2 : Comparison of two giant air showers 


Zenith angle 

Number of water- 
Cerenkov detectors 
and distance range 

Number of 1 m 2 
scintillators and 
distance range 

S (600) m~ 2 

p(600) m~ 2 

p v (600) m 2 

Primary energy: 

Yakutsk calibration 

Hillas relation 


21220296 17683412 

(J Phys G 9 _, 1569 1983) (World Data Catalogue 

pp86-7) 


13 ° 

24 

150 < r < 2170m 
8 

420 < r < 680 m 
157 
64 
66 


35° 

50 

90 < r < 2500 m 


105 

136 


5. 3 x 10 19 eV 

5 . 0 x 10 1 9 eV 1. 1 x 10 2 ° eV 


The estimated error in the assigned size (p(600)) for each of these 
events is ^ 10%; this error includes core location uncertainty, 
stationary error and allowance for lateral distribution uncertainty and 
is so small because of the exceptional symmetry in the detector density 
patterns. The risetime measurements in each event are also in agreement 
with these analyses. Event 17683412 is unquestionably twice as large as 
21220296 which in turn, through the scintillator and water-Cerenkov 
densities, has two independent energy estimates of ^5xl0 19 eV. 


In addition to the 4 events discussed above we have recorded a further 4 
everts which we believe are ^10 2 ° eV. These are not included in our 
energy spectrum because they arrived from zenith angles > 45° and/or the 
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cores fell outside of the array boundary. The flux derived from all 8 


+2 


-1 


x 10 


“16 


m 


1 sr" 1 


events in this total exposure of 657 km 2 sry is 

. V. - ) 

and is consistent with that deduced for the 4 events of Table 1, namely 
I(> 10 20 eV) = [3 *1' 


x 10 1 6 m 2 s 1 sr 1 


4. Discussion and Conclusions . The proven existence of cosmic ray 
events with E > 10*° eV demands explanation. Presumably the source of 
these events must be relatively close to the earth but it can hardly be 
galactic as |b| >40° for all 4 events of Table 1. The inferences drawn 
about the ability of the Cygnus X-3 system to accelerate large fluxes of 
cosmic ray nuclei to 10 1 7 eV/nucleon (Hillas 1984) leads naturally to 
speculation that a suitably scaled up system, perhaps in the nucleus of 
an active galaxy, can accelerate particles to 10 20 eV and beyond. 


Our current best estimates of the integral intensities above 10 18 , 10 19 
and 10 2 ° eV are 


I(> 10 18 eV) = 
I(> 10 19 eV) = 
I(> 10 2 ° eV) = 


(1.9 ± 0.2) 
( 2.1 ± 0 . 2 ) 
(3 ± 2) 


x 10 1 2 m 2 s 1 sr 1 

x 10 1 4 m 2 s 1 sr 1 

x 10 16 m 2 s 1 sr 1 


Further details of our analysis will be published elsewhere. 
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t After preparation of this paper was completed (22 May 1985) the issue 
of Physical Review Letters (22 April 1985), which contains the new 
spectrum deduced from Fly's Eye, reached us. We do not agree with the 
Fly s Eye group's conclusion that an end to the cosmic ray spectrum 
has been observed. A comment on their letter is being prepared for 
submission to Physical Review Letters. 



